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Multi-Phase Flow 



K i 



for each 



Compute 

node of the4-Tet, 

j = 1,2, 3, 4, 5 Where; 



T 



Divide Kj into positive and nega- 
tive values 

Kj 1 " = K J -ifK j i0 f dseO 
Kj = KjifK.<0,elseO 



T 



Algorithm 1. To distribute residuals of 
the finite element operator: (Tj 1 )^^ 



Distribute residual to the nodes 


according to 








k 




^ k J 





Compute tangent operator • 



Assemble system and solve. 



Fig. 8 



